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Grazing scattering of neutral atoms from insulator surfaces is investigated in the intermediate 
velocity range, in which interference effects have been recently observed. To describe this process 
we introduce a distorted-wave method, based on the use of the eikonal wave function, which takes 
into account the phase of the scattering state along the classical projectile path. The eikonal theory 
is applied to evaluate the angular distribution of few keV helium atoms after impinging on a LiF(OOl) 
surface along low-index crystallographic directions. The interest focuses on the role played by the 
projectile polarization produced by cations and anions of the crystal surface. For the considered 
collision system we found a polarization channel, corresponding to the direction <110>, which is 
affected by this effect, while for incidence in the direction <100> the polarization contribution is 
nearly negligible. The proposed eikonal approach, including polarization effects, provides angular 
projectile spectra in good agreement with the experimental data. 

PACS numbers; 



I. INTRODUCTION 

Experimental evidences of interference effects pro- 
duced during grazing scattering of fast atoms from insula- 
tor surfaces were recently presented in Refs. [1-3]. Under 
axial surface channeling conditions, reported measure- 
ments of scattered projectile distributions display well 
defined spots, associated with diffraction patterns orig- 
inated by the periodic structure of the crystal. Even 
though the diffraction of particles from crystal surfaces 
has been a well-understood phenomena from the begin- 
ning of quantum mechanics, the importance of these ex- 
perimental results is due to the effect not being expected 
to be observable for light atoms with energies in the keV 
range, whose de Broglic wavelengths arc some orders of 
magnitude smaller than the shortest interatomic distance 
in the crystal. 

Two different interference mechanisms have been pro- 
posed to explain the experimental observations. The first 
of them [1, 2] is related to the diffraction by a periodic 
lattice and relies on the assumption that the projectile 
motion perpendicular to the axial incidence channel can 
be decoupled from the parallel one. The semiclassical ap- 
proach [2] proposed for the description of this mechanism 
predicts a maximum in the distribution when the com- 
ponent normal to the scattering plane of the final projec- 
tile momentum coincides with a reciprocal lattice vector. 
Such a result, also confirmed by means of a wave packet 
propagation approach [2], is in accord with the experi- 
mental data for the lowest impact energies [1, 2]. The 
second mechanism [3], called supernumerary rainbow, is 
originated by the corrugation of the surface potential, 
which gives rise to a quantum interference between pro- 
jectiles emerging from the surface with the same direction 
but reflected at different turning points [4] . 

To describe the experimentally observed patterns, in 



this article we introduce a distorted-wave model, named 
surface eikonal approximation, which is valid for small 
de Broglic wavelengths of incident atoms. This method 
makes use of the eikonal wave function [5] to represent 
the elastic collision with the surface, while the projec- 
tile movement is classically described, taking into ac- 
count different initial conditions. The surface eikonal 
approach can be considered as an extension of the well- 
known Glauber approximation [6] for collisions with cor- 
rugated surfaces [7] instead of atoms, but considering ax- 
ial channeled trajectories. It includes both interference 
mechanisms and the idea behind it essentially coincides 
with that of the semiclassical formalism [4] used in Ref. 
[3]. 

The surface eikonal approach is here employed to de- 
scribe angular distributions of swift He" atoms scattered 
off from a LiF(OOl) surface, for which there are experi- 
mental data available [1-3]. As the considered process is 
very sensitive to the description of the surface potential, 
the aim of the work is to investigate the influence of the 
polarization on the interference patterns. In our model 
the interaction of the incident atom with the crystal sur- 
face is represented as a sum of individual interatomic 
potentials, which take into account the contribution of 
the different ionic centres of the insulator material [8]. 
To evaluate the interatomic potentials we use the Abra- 
hamson approximation [9] , adding the asymptotic contri- 
bution of the projectile polarization. The role of the po- 
larization is analyzed for incidence along the <100> and 
<110> channels, finding that polarization effects are im- 
portant for this latter crystallographic direction. Atomic 
units (e^ = h = irif. = 1) are used unless otherwise stated. 
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II. SURFACE EIKONAL APPROXIMATION 

Let us consider the grazing impact of an atomic pro- 
jectile (P), with mass mp, on a crystal surface (S). As a 
result of the collision, the projectile with initial momen- 
tum Ki is elastically scattered from the surface, ending 
in a final state with momentum K f. The frame of refer- 
ence is fixed on a target ion belonging to the first atomic 
layer, with the surface contained in the x — y plane and 
the z versor perpendicular to the surface, aiming towards 
the vacuum region. 

We assume that the state associated with the col- 
lision system satisfies the time-independent Schrodinger 
equation for the Hamiltonian 



H = — 



1 



2mp 



Vsp{Rp), 



(1) 



where Rp denotes the position of the center of mass of 
the incident atom and Vsp is the surface-projectile in- 
teraction. As initial condition, when the projectile is 
far from the surface, "iff tends to the state with 
<^>j{Rp) = (27r)-3/2exp(iJ?j • Rp), j = i{f), the initial 
(final) unperturbed wave function. 

The central magnitude to describe the elastic scatter- 
ing process is the transition matrix, which reads 



(2) 



T,f = / dRp ^f{Rp) VspiRpMiRp)- 



In the energy range of interest, Eq. (2) can be expressed 
in terms of the classical trajectory of the projectile - TZp 
- by means of the substitution Rp = TZp, like in the 
usual semiclassical formalism [10]. The position TZp of 
the incident atom at a given time t is governed by the 
Newton equations associated with the potential Vsp, ver- 
ifying the relation 'JZp{Ros,t) = Ros + ZqZ + j*i^dt' 
v{Ros,t'), where v{Ros,t) is the classical velocity of the 
projectile, Ros = (^0,^0,0) identifies its initial position 
on the surface plane and Zo — > -l-oo. A sketch picture 
of the projectile path and the coordinate system is dis- 
played in Fig. 1. By replacing the integration variables 
Rp = {Xp,Yp, Zp) by the new ones {Xo,Yo,t} in Eq. 
(2), the transition matrix is expressed as [11] 



+ 00 

dRos I dt 



V^{Ros,t) 



(3) 



where Vz{RosTi) is the component of v normal to the 
surface. 

Since the de Broglic wavelength of the incident pro- 
jectile, A = 2-K/Ki, is sufficiently short compared with 
the characteristic distance of the surface potential, we 



approximate the scattering state by means of the 
eikonal wave function [5], i.e. 

*+(7?p) ~ xr^iUp) = $,(7?p)exp(-z,7(7?p)), (4) 

where ?7(7?p) is the eikonal phase, defined as 

t 

Tj{np{Ros,t)) = [ dt' Vsp{np{Ros.t')). (5) 



By introducing the function Xi in Eq. (3) the eikonal 
transition matrix reads 



+ 00 



(eik) 



dRos / dt 



V:,{Ros,t) 



(27r)3 

— OO 

exp[-zQ • .7?p - z?7(7?p)] Vsp{np), (6) 

where Q = Kf — Ki is the projectile momentum transfer 
and the final momentum Kf satisfies the energy conserva- 
tion, i.e. Kf = Ki. The differential probability, per unit 
of surface area, of elastic scattering with final momentum 
Kf in the direction of the solid angle / is obtained from 



{eik) 
if 



where Tj^p''^ dc- 



Eq. (6) as dP/d^f = {2Tr)*ml 

notes the eikonal T-matrix element, normalized per unit 
area. Note that the main difference between the usual 
eikonal scattering amplitude [5, 7] and Eq. (6) arises 
from the use of axial channeled trajectories instead of 
straight-line ones. 

The first Born T-matrix element can be derived from 
Eq. (6) by neglecting the eikonal phase; that is, by fixing 
?7(7?p) = 0. It reads 



{Born) 



1 



+ 00 



dRo 



dt 



VziRos,t) 



(27r)3 

— OO 

exp(-iQ • 7?.p) Vspilip)- 
III. INTERACTION POTENTIALS 



(7) 



In this work, the projectile-surface potential Vsp con- 
tains the static and polarization interactions, i.e., Vsp 



(st) 
SP 



y^p°'\ Due to the insulator character of the sur- 



face, both term can be derived by considering the surface 
as composed by independent target ions. Consequently, 
the static potential Vgp^ is expressed as a sum of in- 
dividual interatomic potentials, Vst{R), which represent 
the static interaction of the incident atom with solid ions 
placed at different lattice sites [8] . Following the Lenz en- 
ergy functional [12, 13], for the two types of target ions 
- alkali-metal and halide- the static ion-atom interaction 
is found to be the sum of three terms, 



Vst{R) = VcouiiR) + V^niR) + V^,h{R). (8) 
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The first term is the well-known electrostatic Coulomb 
interaction, 

VcouiiR) = \j j drdT^Drir^y^^DpiT^ - R), (9) 

where Dxir) = ZT5{r)~ and Dp{r) = Zp5{r) — 

Pp{r) are the target and projectile charge densities, 5 is 
the Dirac delta fmiction, situated at the position of the 
nucleus, Zt (Zp) is the target (projectile) nucleus charge, 
and prp (pp) is the target (projectile) electronic density. 
Note that Vcoui is composed by four terms, including 
the internuclear and electron-electron repulsions as well 
as the attractive electron-nucleus potentials. 

By employing the Abrahamson approximation [9] the 
second term of Eq.(8), named the kinetic potential, reads 



of the radius of the atom, that is, i?oi = ('')Ti + ('')p! 
where (f)^- {{r)p) is the target (projectile) mean ra- 
dius. We employed the values 1.09, 0.67, and 1.41 a.u. 
for the He, Li+ and F~ mean radii, respectively. Far 
from the surface, the diagonal terms (i = j) of the po- 
larization potential given by Eq. (12) satisfy the well 
known behavior ~a/2R'^, but the extra-diagonal terms 
{i ^ j) arc weighted by a directional factor that depends 
on the crystal ordering [15, 16]. Notice that as we are 
dealing with neutral projectiles, we have not taken into 
account the dynamic polarization of the surface ions [8] 
because this effect represents a higher-order correction of 
the interatomic potential (cc R^^). 



IV. RESULTS 



2.871 



= / dr 



Prir) + Pp{r-R) 



drp'^''{r) - / df pTi^-R) 



5/3 



5/3, 



(10) 



This potential is essentially positive and represents the 
reaction to the compression of the electronic density, con- 
sidered as a free-electron gas. The third term describes 
the exchange potential within the local approximation 
and it reads 



V^ch{R) 



-1) 



0.738 



dr 



Prir) + Pp{r-R) 



4/3 



(11) 



df py,^{r) — / df p'^p'^{f—R) 



4/3^ 



Finally, the potential Vjp"'^ takes into account the polar- 
ization of the neutral projectile in the presence of target 
ions, which is not included in the original Abrahamson 
model. Following the usual derivation of the atomic po- 
larization potential [5], although in this case, not for only 
one perturbative charge but a collection of them, which 
represent the different target ions, the asymptotic polar- 
ization potential reads 



v^'/\r) 



7(00) 



2^[Rl, + R}) 



{R,-R,) 



7(00) 



(Rf 



R]y 



(12) 

where the sum formally includes all the target ions of 
the crystal, a is the polarizability of the projectile, with 
a = 1.38 a.u. for Helium [14], and Ri (Rj) represents the 
position vector of the projectile with respect to the target 
ion labelled as i (j), with R^ = R^/R^. In Eq.(12), Z^^ 
is the residual charge of the target ion at long distances, 
being z!^^ = 1 for Li+ and z!f^^ = -1 for F^. At short 
distances the polarization contribution of the target ion 
i is reduced with a cutoff, which is always of the order 



We applied the model to neutral helium atoms imping- 
ing grazingly on a LiF crystal surface under axial surface 
channeling conditions. The impact energy ranged from 
0.2 to 8.6 keV, corresponding to the experiments of Refs. 
[1-3]. In the crystal surface, ions belonging to the top- 
most atomic layer were slightly displaced from their equi- 
librium positions, in accord with Ref. [17]. 

To describe the projectile-surface potential we em- 
ployed the punctual model of Ref. [8], evaluating the 
He-Li+ and He-F~ interatomic potentials from Eq.(8). 
Hartree-Fock Slater wave functions from Clementi-Roetti 
[18] were used to calculate the electronic densities 
and Pp . It allowed us to derive a closed form for Vcoui , 
while Vkin and Vxch were obtained from numerical in- 
tegrations. In Fig. 2 it seemed convenient to plot the 
scaled expression W{R) = V{R) * R{1 + 2R^) for Li+ 
and F~, respectively, where V{R) includes the static po- 
tential [Eq.(8)] plus the diagonal polarization contribu- 
tion, i.e. the i = j term of Eq.(12). From the figure we 
can differentiate two different regions of the interatomic 
potentials. As i? ^ 0, W{R) ZtZp and the sharp in- 
crease at the origin corresponds to the electrostatic con- 
tribution VcouiiR)i while the maximum at intermediate 
distances is mainly due to the statistical contribution, 
i-e. Vkin{R) + Vxch{R)- Note that present static poten- 
tials are almost indistinguishable from the ones of Gor- 
don and Kim [19] (empty circles), employed in Ref. [3]. 
The asymptotic limit of as i? 00 is affected by the 
polarization, i.e. V * R{\ + 2i?^) —a. 

The projectile trajectory was derived from classical 
dynamics with the Runge-Kutta method. At every step 
we took into account the 4*'' order nearest neighbor tar- 
get ions (i.e. 8x8x4), which includes the interaction 
of the projectile with the topmost atomic layer and three 
more layers below it. We have made sure our results do 
no depend on the considered number of nearest neigh- 
bors by increasing this number to include up to 8*'' order 
nearest neighbors (i.e., 8 atomic planes). 

The evaluation of the eikonal transition matrix involves 
an integration on the starting point i?os of the classical 
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trajectory, which was calculated with the MonteCarlo 
technique, varying Ros on the area of the unit cell as 
a consequence of the surface invariance. In every case we 
considered around 10^ classical trajectories with random 
initial positions, and this number was varied in order to 
test the convergency of our calculations. The further in- 
tegration on t involved in Eq. (6) was numerically solved 
with a relative error lower than 0.1%. To obtain the 
differential probability dP/dflf, we have to add the T- 
matrix elements corresponding to different values of Ros 
that lead to the same final momentum Kj. For this pur- 
pose we employed a grid for the angles 0f and of 
100 X 100 points, where 9f and (p^ are the final polar and 
azimuthal angles, respectively, of the final momentum 
Kf. In all calculations we oriented the x versor along 
the low-index direction of the crystal surface coinciding 
with the impact direction; therefore, the azimuthal angle 
ipj: is measured with respect to the incidence direction 
on the surface plane (see Fig. 1). In full accord with 
the experiments of Rcfs. [1-3] we found that under axial 
surface channeling conditions the relation O'j + ip'^j: ~ 9^ is 
almost strictly verified by all classically scattered projec- 
tiles and consequently, the angular projectile distribution 
shows the usual banana shape [20] . 

We start the analysis by considering the experimental 
case of Fig. 5 of Ref. [1]; that is, 3 keV ^He atoms im- 
pinging on a LiF(OOl) surface along the crystallographic 
direction <110> with a glancing angle {di = 1.1 deg). 
This collision system looks adequate for the eikonal de- 
scription because the de Broglie wavelength of the in- 
cident atom (A = 0.0057 a.u.) is almost three orders 
of magnitude smaller than the characteristic interatomic 
distance. In Fig. 3 we plot the differential probability 
dP/dipj:, as a function of the azimuthal angle ipj, mul- 
tiplying the results by an arbitrary factor in order to 
show the different curves separately. The eikonal spec- 
trum displays strong interference signatures, presenting 
pronounced maxima symmetrically placed with respect 
to the incidence direction, which corresponds to 1^9^ = 0. 
This interference pattern can be directly compared with 
the experimental spots of Ref. [1], which are displayed 
with stars, numbering them from the central one. The 
eikonal distribution nearly agrees with the experimen- 
tal one, although the eikonal maxima associated to the 
peaks ±1, ±2, ±3, and ±4 are slightly shifted to higher 
values. Notice that the extreme angles of the eikonal 
spectrum are related to the rainbow scattering and the 
corresponding maxima display a sharp shape. These 
peaks are also present in the classical distribution, de- 
fined as the number of projectile trajectories reaching a 
given final azimuthal angle which is shown in an abso- 
lute scale in Fig. 3. The classical scattering distribution 
presents the typical rainbow profile [21], with only two 
maxima around the extremes of the angular spectrum - 
the rainbow angles. Then, the absence of intermediate 
structures in the classical spectrum confirms the concept 
that interference effects are a consequence of quantum co- 
herence between projectiles moving along different paths 



but ending in the same final state. 

With the aim of analyzing the influence of the polariza- 
tion of helium atoms, in Fig. 3 we also plot eikonal values 
obtained by neglecting the polarization potential; that is, 
by dropping Vg^p^^ in the projectile-surface interaction. 
We found that for incidence along the direction <110>, 
the angular distribution of scattered atoms is affected 
by the projectile polarization. When the polarization is 
not included in the calculation, the central maximum be- 
comes a minimum, modifying the total number of peaks 
displayed by the eikonal distribution. In turn, the ex- 
treme maxima, associated with the rainbow angles, are 
only slightly altered by the polarization. Both angular 
regions - central and external- of the eikonal spectrum 
are associated with different zones of the interatomic po- 
tentials that are probed by axial channeled projectiles. 
He" atoms that reach azimuthal angles (pf near move 
over the ionic rows that form the channel, farther than 2 
a.u. from the surface, interacting with the long-distance 
contribution of the surface-projectile potential. As a such 
contribution is dominated by the term corresponding to 
the polarization potential, given by Eq.(12), it explains 
the inffuence of this effect on the central zone of the spec- 
trum. Projectiles that end in the rainbow angular region, 
instead, suffer closer collisions with F ionic centers, being 
affected by the short-distance behavior of the interatomic 
potentials, which is determined by coulombic and statis- 
tical contributions. 

In Fig. 4 we investigate the elastic scattering along the 
direction <100> by considering a higher impact energy 
(8.6 kcV). The eikonal differential probability is plotted 
in Fig. 4 (a) as a function of the final azimuthal angle, 
together with experimental spots of Fig. 2 of Ref. [3]. 
For this collision system, in addition to the two rainbow 
maxima, the eikonal distribution presents four similar 
peaks, symmetrically placed around ipj = 0, and a very 
small central maximum. The number of main maxima 
of the eikonal profile coincides with that of the experi- 
mental pattern [3], although the positions of the peaks 
are again shifted to higher values in comparison with the 
experimental ones. In all the cases we found that slight 
changes in the interatomic potentials produce substantial 
modifications in the angular spectrum of scattered pro- 
jectiles. Hence, discrepancies between theoretical and ex- 
perimental spectra could be associated with very subtle 
differences in the projectile-surface potential. In Fig. 4 
(b) we compare experimental intensities [3] with eikonal 
probabilities, now plotted in linear scale, as function of 
the deflection angle O, dcflned as = arctan((y9j /Of). 
Taking into account that our theoretical results were ob- 
tained by considering fixed positions of the target ions, 
without including the thermal vibration, and they were 
not convoluted with experimental conditions, the eikonal 
model reproduces fairly well the main features of the ex- 
perimental spectrum. 

In order to investigate the effect of the polarization in 
the channel <100>, in Fig 4 (a) we show eikonal values 
derived by eliminating the polarization potential. Re- 



5 



markably, cikonal results with and without including the 
projectile polarization agree with each other for incidence 
along the <100> direction, indicating that polarization 
effects play a minor role in this channel. It is a conse- 
quence of the ordering of the halide and alkali ions in- 
volved in the axial surface channeling. As observed from 
Eq. (12), when the projectile moves along the channel 
far from the surface plane, the factors of the polariza- 
tion potential coming from F~ and Li"*" have opposite 
signs and they compensate their contributions to order 
when F~ and Li+ ions are placed in front of each 
other, as it happens in the <100> direction. Further- 
more, within a row model, the <100> rows - formed 
by alternate cations an anions - display a neutral charge, 
which reduces the polarization of the incident atom. In 
the <110> direction, instead, not only arc there sepa- 
rated cation and anion rows, with positive and negative 
net charges respectively, but also Li+ and F^ions are not 
in front of each other along the channel, which originates 
an effective polarization potential. This is the reason why 
polarization effects become evident for incidence along 
the <110> direction but not in the channel <100>. 

Besides, in Fig. 4 (a) we also show the angular dis- 
tribution obtained within the first Born approximation 
[Eq.(7)], which is derived from Eq. (6) by eliminating 
the eikonal phase. The Born profile displays a different 
diffraction pattern, with a broad central maximum, not 
present in the experiment, indicating that interference 
structures of the surface eikonal model are affected by 
the phase rj, given in Eq. (5). However, note that differ- 
ences between cikonal and Born distributions vary with 
the considered collision system. 

Finally, in Fig. 5 we considered the incidence condi- 
tions of Ref. [2] , which correspond to a smaller impact 
energy (0.2 keV). Notice that this energy is close to the 
limit of validity of the eikonal model, which is expected 
to be adequate for high velocities. For scattering along 
the direction <110>, the eikonal differential probability 
is displayed in Fig. 5, as a function of the azimuthal an- 
gle, comparing it with the spots of Fig. 1 of Rcf. [2]. 
Also in this case, the agreement of the cikonal theory 
with the experiment is reasonable good. Both profiles 
- eikonal and experimental - present similar structures, 
with a central maximum and two additional peaks, not 
equally spaced, to each side. However, the experimen- 
tal peaks ±1 are narrower than the eikonal ones, and 
small structures around the rainbow angles are absent 
in the theory, corresponding to the worst disagreement 
found in the present work. Discrepancies between the 
theory and the experiment can be again attributed to 
extremely subtle distinctions in the projectile-surface po- 
tential. Moreover, we should mention that the small rum- 
pling {d = 0.037 a.u.) of the surface ions introduced in 
our model [17] affects the interference pattern. 

Again, like in Fig. 3, the central zone of the eikonal 
spectrum of Fig. 5 is associated with the long-distance 
behavior of the surface interaction, which is governed by 
the projectile polarization. When Vg^''^ is dropped, the 



central maximum of the cikonal distribution completely 
disappears, in disagreement with the experimental data. 

To investigate in detail the central zone of the eikonal 
spectrum, in Fig 6 we plot the first ten projectile trajec- 
tories, provided by the MontcCarlo code, that contribute 
to the distribution at the final azimuthal angle (pj: = 0. 
For the collision system of Fig. 3 we observe that all the 
atoms that end in this angular region move just over F~ 
or Li^ rows. In this case, turning points corresponding 
to the z- movement are almost independent of the mo- 
tion perpendicular to the scattering plane, being approx- 
imately situated 2.8 a.u. (2.2 a.u.) above the topmost 
atomic layer for projectiles moving over F~ (Li+) rows. 
From Fig. 6 (c), the transversal kinetic energy, defined 
as ii;^'™^ = mp{vy + ^1)12, slightly increases just before 
and after reaching the collision region, indicating that 
incident atoms are affected by an attractive polarization 
potential. The projectile-surface potential along classical 
trajectories, shown in Fig. 6 (d), displays an oscillatory 
pattern produced by the interaction with the different 
ionic centers of the crystal surface. Consequently, the to- 
tal transversal energy E_\_ = + Vsp presents fluctu- 
ations along the classical projectile path. However, the 

mean value {E±^ = ^iJ^*^™'^ + {Vsp) keeps equal to 

the initial value Ei^ = mpvf^/2 along the whole trajec- 
tory, supporting to some extend the decoupling of the 
transversal movement from the parallel one, proposed in 
Ref. [2]. 



V. SUMMARY 

In conclusion, we have developed a surface eikonal ap- 
proach to deal with interference patterns produced by im- 
pact of swift atoms on insulator surfaces. The proposed 
method has been applied to few keV He atoms grazing 
impinging on LiF(OOl) along the <110> and <100> di- 
rections. Projectile spectra derived with the eikonal ap- 
proximation display well-defined interference structures, 
originated by atoms that follow different paths but end 
scattered with the same final momentum. As the pro- 
jectile distribution strongly depends on the description 
employed to represent the projectile-surface interaction, 
the study focused on the influence of the projectile polar- 
ization on the angular spectrum. We conclude that the 
polarization potential is essential to describe the elastic 
scattering along the <110> channel, while in the direc- 
tion <100> its contribution is negligible. Angular spec- 
tra derived from the cikonal model, including the po- 
larization effect, are in concordance with the available 
experimental data [1-3] . But a better representation of 
the surface potential, taking into account that target ions 
are part of a surface, might modify the present results. 
Then, by including more precise electronic densities this 
method may be useful to investigate very delicate items, 
such as long-distance potentials or crystal ion displace- 
ments, which are difficult to make evident experimentally 
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FIG. 1: Schematic depiction of the coordinate system. 



FIG. 2: Scaled interatomic potentials for (a) He-Li+ and (b) 
He-F~. Solid line, potential including the diagonal (i—j) po- 
larization contribution; dashed line, static potential, without 
polarization; and empty circles, results reported by Gordon 
and Kim in Ref. [19]. 



FIG. 3: Azimuthal angular distribution of elastic scattered 
projectiles for 3 keV ^He atoms impinging on LiF(OOl) along 
the direction <110>, with the incidence angle 6i — 1.1 deg. 
Solid line, differential probability derived from the surface 
eikonal approach, including polarization effects; dashed line, 
surface eikonal results without including the projectile polar- 
ization. Full stars, experimental spots of Fig. 5 of Ref. [1]. 
Empty circles, classical distribution, as explained in the text, 
in absolute scale. 



FIG. 4: Similar to Fig. 3 for 8.6 keV *He atoms impinging on 
LiF(OOl) along the direction <100>, with the incidence angle 
9i = 0.71 deg. Dash-dotted line, first Born approximation 
[Eq. (7)]. (a) Full stars, experimental spots, and (b) thick 
solid line, experimental intensity, both drawn from Fig. 2 of 
Ref. [3]. 



FIG. 5: Similar to Fig. 3 for 0.2 keV *He atoms impinging 
on LiF(OOl) along the direction <110>, with the incidence 
angle 6i = 1.5 deg. Full stars, experimental data from Fig. 1 
of Ref. [2] 



FIG. 6: Classical projectile trajectories ending with a final 
azimuthal angle ip^ near to 0, as a function of the coordi- 
natate Xp parallel to the axial channel, for the collision sys- 
tem of Fig. 3. (a) Coordinate Yp perpendicular to the scat- 
tering plane; (b) coordinate Zp perpendicular to the surface; 
(c) transversal kinetic energy along the trajectory, defined as 
E± = mp{Vy + v1)/2\ (d) projectile-surface potential along 
the trajectory. 
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